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Two candidates for "almost-invariant" toroidal surfaces passing through magnetic islands, namely quadratic-flux-minimizing 
(QFMin) surfaces and ghost surfaces, use famiUes of periodic pseudo-orbits (i.e. paths for which the action is not exactly extremal). 
QFMin pseudo-orbits, which are coordinate-dependent, are field lines obtained from a modified magnetic field, and ghost-surface 
pseudo-orbits are obtained by displacing closed field lines in the direction of steepest descent of magnetic action, ^ A-dl. A gener- 
alized Hamiltonian definition of ghost surfaces is given and specialized to the usual Lagrangian definition. A modified Hamilton's 
Principle is introduced that allows the use of Lagrangian integration for calculation of the QFMin pseudo-orbits. Numerical calcula- 
tions show QFMin and Lagrangian ghost surfaces give very similar results for a chaotic magnetic field perturbed from an integrable 
case, and this is explained using a perturbative construction of an auxiliary poloidal angle for which QFMin and Lagrangian ghost 
surfaces are the same up to second order. While presented in the context of 3-dimensional magnetic field line systems, the concepts 
are applicable to defining almost-invariant tori in other 1 ^ degree-of-freedom nonintegrable Lagrangian/Hamiltonian systems. 
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1. introduction 

The understanding of nonintegrable Hamiltonian systems is 
greatly simplified if one can construct a coordinate framework 

L^Zj^based on a set of surfaces that are either invariant under the 
dynamics or, where this is impossible, surfaces that are almost- 
invariant. As invariant tori and cantori in nonintegrable sys- 
tems can be approximated by sequences of periodic orbits, the 
theory of almost-invariant surfaces is built around periodic or- 
bits, which constitute the remanent invariant sets surviving af- 
ter integrability is destroyed by symmetry-breaking perturba- 
tions. We consider two classes of almost-invariant surfaces, 
quadratic-flux-minimizing (QFMin) surfaces [1] and ghost sur- 

0^ faces [2, 3]. 

^P, Almost-invariant tori are important in the theory of magnetic 
>• confinement of toroidal plasmas, in particular to the theory of 

*kp transport in chaotic magnetic fields [4], and we set this paper in 
the context of the nonintegrable magnetic fields, B, encountered 
in devices without a continuous symmetry. However, as mag- 
netic field lines are orbits of a 1 ^ degree-of-freedom Hamilto- 
nian system, [5] the discussion is applicable, with appropriate 
translations of terminology, to any such system — e.g. in this 
paper we use "magnetic field line" and "orbit" interchangeably. 

In Sec. 2 we introduce our general, arbitrary background 
toroidal coordinate system s, 0, and an auxiliary poloidal an- 




Figure 1 : A sketch of the general curvihnear toroidal coordinate system de- 
scribed in the text. 
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gle &(s, 6, that allows us to define the quadratic flux in a form 
independent of the choice of 6. In Sec. 3 we introduce the mag- 
netic action integral, its first and second variations and Hamil- 
ton's principle, while in Sec. 4 we introduce QFMin and (gen- 
erahzed) ghost-surface pseudo-orbits as alternative strategies 
for continuously deforming the action-minimax orbit associated 
with an island chain into the corresponding action-minimizing 
orbit. 

In Sec. 5 we present numerical results for field-line Hamilto- 
nians of the form;t'o(iA) + ^Xi ('A' 0^ where the flux function i// 
plays the role of a momentum canonically conjugate to 0, and 
e parametrizes the strength of the perturbation away from the 
integrable case described by the action-angle Hamiltonian xo- 
Plots are presented comparing the uncorrected (i.e. with @ = 6) 
QFMin and Lagrangian ghost curves of Ref. 2,superposed on 
field-line puncture plots in a Poincare surface of section. Two 
cases with dififerent strengths of perturbation are shown, both 
quite strongly chaotic and both showing that the differences be- 
tween even uncorrected QFMin and ghost curves are very small 
(except for some higher-order surfaces, in the more strongly 
chaotic case). This suggests that the two, seemingly very differ- 
ent, approaches to defining almost-invariant tori may be unified 
by appropriate choice of 0, and that this will differ from 6 by 
an amount small in e. 

In Sec. 6 we introduce a modified form of Hamilton's Prin- 
ciple that gives QFMin pseudo-orbits as extremizers of a pseu- 
doaction. Section 7 gives the canonical, Hamiltonian form of 
this action principle, while Sec. 8 discusses the transformation 
to the Lagrangian form. In Sec. 9 we derive a consistency con- 
dition that must satisfy for corrected QFMin surfaces to be 
Lagrangian ghost surfaces, finding in Sec. 10 an expression for 
a choice of the auxiliary angle that satisfies this criterion up 
to first order in e. The difference between uncorrected QFMin 
and ghost/corrected-QFMin pseudo-orbits is shown indeed to 
be very small, O(e^). 

In Sec. 1 1 we sketch our finite-element variational method 
for numerical construction of QFMin surfaces using the new 
Hamilton's Principle introduced in Sec. 6, and in Sec. 12 we dis- 
cuss the numerical construction of ghost surfaces via Galerkin 
projection onto the finite element basis. 

Appendix A contains a derivation of the Euler-Lagrange 
equation for QFMin pseudo-orbits in the canonical represen- 
tation, and Appendix B shows the relation between the gener- 
aUzed definition of ghost pseudo-orbit given in Sec. 4 and our 
more standard Lagrangian form [2], used in the numerical work 
and in Sec. 9. 

2. Coordinates and fluxes 

As depicted in Fig. 1, we assume a general, essentially arbi- 
trary curviUnear toroidal coordinate system s(r), 6(r), ^(r) has 
been established, where r is a point in Euclidean 3-space and 
6 and f are respectively poloidal and toroidal angles label- 
ing points on the toroidal isosurfaces of s, nested around the 
curve along which 6 is singular (s increasing outward). We as- 
sume the nonorthogonal basis {e', e*, e^ ) = {Vs, V6, V^) is right 
handed, as is its reciprocal basis {e^, e^, e^ } = {5,r, dot, d^r). 



The directed infinitesimal area element on an arbitrary sur- 
face F is JS = dOd^njn'VdxV^, where n is the unit normal 
at any point on F. Thus the net magnetic flux crossing an ar- 
bitrary torus F (which we assume to contain the 0-coordinate 
singularity curve) is 



n-B 



(1) 



n-V6»xVf 

This integral is independent of choice of coordinates. In fact 
the absence of magnetic monopoles impUes that ipi vanishes 
identically, so it is independent of the choice of F also, whether 
it be a magnetic surface (invariant torus of the field-line flow) 
or otherwise. 

Thus, to measure the amount by which F departs from being 
a magnetic surface, we are led to define the positive definite 
quadratic flux [1], defined with the aid of a new generalized 
poloidal angle @{s, 6, f ), 
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The quadratic flux ^2 is independent of the choice of base co- 
ordinates s, 0, f , but depends on the choice of 0. 

In the numerical work presented in this paper, has been 
chosen equal to the given angle 0. However in the formal de- 
velopment we distinguish it from so we can explore the con- 
sequences of making different choices, in particular whether it 
can be chosen so that QFMin tori coincide with ghost tori. 

3. Magnetic action integral 

The field-Une action S [6] is a functional of a path C in Eu- 
clidean 3-space, points on which we take to be labeled by the 
toroidal angle ^, which thus takes on the role played by time 
in a more conventional Hamiltonian system. In this paper we 
confine our attention to paths that are closed loops, with in- 
creasing by 2np when f increases by 2nq (p and q > being 
mutually prime integers), so the average rate of increase of 
along the path is the rational fraction p/q, where the angular 
frequency t is called the rotational transform. 

The magnetic action is defined by 

S[C]= I A-dl= I A-rd^, (3) 
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where the single-valued function A(r) is a magnetic vector po- 
tential for the magnetic field, B = VxA, and dl = td( is an 
infinitesimal Une element tangential to C. A superscript dot de- 
notes the total derivative with respect to ^, so that r-Vi^ = 1. 
Hamilton's Principle is the statement that S is stationary, with 
respect to variations 5r of C, when C is a segment of a physical 
orbit (in our case a magnetic field line). If C is an open segment 
the variations are to be taken holding the endpoints fixed, but if 
(as we assume) C is a closed loop then the variations are uncon- 
strained because the endpoint contributions cancel. Then, after 
integration by parts, we have the expansion for the total change 
in S 
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where the first functional derivative is given by 

rxB . 
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Hamilton's Principle is now readily verified: The Euler- 
Lagrange equation SS/6r = is satisfied if r = B/B^, i.e. on a 
magnetic field line. 

The symmetrized Hessian operator is 



I = Bxl - IxB — 

6r6r 

+ rx(VBf - (VB)xr , 



(6) 



where I = e^e ' + Cge* + = e^e^ + e^e^ + e^Cf is the identity 
dyadic and superscript ^ denotes the transpose. 

Note also that variations 6r - r(^ + 6^ - r(f) = r + r6(, + 
^r(5^)^ + . . . that simply relabel the path can be verified to leave 
S invariant for arbitrary 6^{^), as expected. Thus, to find a 
unique minimizer, we must suppress this degree of freedom in 
the allowed variations. To this end we constrain 6r to the tan- 
gent plane of the poloidal surface of section f = const at r, 
denoting the constrained variation by Ar = ts6s + teSO. (Pro- 
vided B, Cj and eg are linearly independent, the two compo- 
nents of the Euler-Lagrange equation, 6S/Ss = e^-rxB = 
and 6S/66 = e^-rxB - imply rxB = 0, so the third compo- 
nent of the Euler-Lagrange equation, 6S/6( = e^-rxB = 0, is 
redundant.) 



4. Ghost and QFMin surfaces 

In an integrable system a continuous family of (p, g')-periodic 
orbits, each extremizing the action, exists, defining an invari- 
ant torus with rotational transform p/q. However such invari- 
ant tori are not structurally stable — small perturbations to the 
system destroy integrability (see e.g. Fig. 2), leaving only iso- 
lated action-extremizing periodic orbits. In fact (assuming a 
twist condition holds), by the Poincare-Birkhoff theorem [7], 
only two distinct periodic orbits survive in a given (p, q) island 
chain, namely the action-minimizing orbit [8], and an action- 
minimax orbit. The minimizing orbit is a hyperbolically un- 
stable "X-point" orbit in the chaotic separatrix region of the 
island chain while the minimax orbit threads the centers of the 
islands. The minimax orbit may be elliptically stable, or, after 
a period-doubling bifurcation, become unstable but continue as 
a {p, q) hyperbolic orbit accompanied by a daughter {2p, 2q) el- 
liptic orbit that does not affect the following discussion. An 
almost-invariant torus, the "ghostly remnant" of an invariant 
torus, is formed from a family of pseudo-orbits, labeled by a 
parameter t, that "fill in" the regions between the minimizing 
and minimax orbits so as to form a torus F. 

The pseudo-orbits of a ghost torus are defined [2] by deform- 
ing the minimax orbit via an action-gradient flow 
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where D/Dt denotes the total r-derivative at fixed i.e. 
D/Dt = {Ds/DT)ds + (DO / DT)dg, and Pghost is a symmetric 
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Figure 2: A comparison of the intersections of unconnected QFMin surfaces 
(tliick daslied lines) and Lagrangian gliost surfaces (thin lines) with the sur- 
face f = 0, with the choices s = ip, & = 8 described in Sec. 7. The two 
almost-invariant surface definitions are almost indistinguishable in this moder- 
ately chaotic case, described in Sec. 5. A Poincare plot (red dots) is also shown. 



nonnegative dyadic. The most natural choice of Pghost might 
seem to be Ppoi s I - e^e^/|e^p projecting onto the poloidal tan- 
gent plane at r, as this is independent of the choice of s and 6. 
However, as we show in Appendix B, this does not correspond 
with that required to recover the usual Lagrangian definition [2] 
of ghost surfaces. 

At the periodic orbits, the action gradient is zero. Begin- 
ning from the minimax orbit, we initially push the curve in the 
decreasing direction (provided by the eigenfunction of the Hes- 
sian with negative eigenvalue) and then evolve the curve ac- 
cording to the action gradient flow Eq. (7). 

Using Eq. (7) in Eq. (4) we find 
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so the sequence of pseudo-orbits tends monotonically toward 
the minimizing orbit, tracing out a surface, which we call the 
ghost surface. 

Ghost surfaces display several attractive properties (proved, 
for Lagrangian ghost orbits, in the case of symplectic maps [3]; 
verified numerically for continuous time/magnetic field systems 
[2]). Their intersections with a surface of section are graphs 
when plotted in canonical phase-space coordinates, so that each 
line of constant crosses the ghost curve only once. They are 
nonintersecting, and thus a discrete selection of ghost surfaces 
may be used as a framework for a generalized action-angle-like 
coordinate system for chaotic fields. Furthermore, there ap- 
pears to be a close correspondence between ghost surfaces and 
the isotherms resulting from strongly anisotropic heat-transport 
in chaotic magnetic fields [4]. However, as it stands, there is a 
significant disadvantage to their construction: the gradient flow 
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Figure 3: A comparison of unconnected QFMin curves (tliick lines) and ghost 
curves (thin lines) for a more strongly chaotic case described in the text. Some 
cases where QFMin curves violate the graph property are seen. 

vanishes as one approaches integrability. In this limit, the con- 
struction of the ghost surfaces becomes arbitrarily slow! 

A QFMin surface is one that minimizes (p2 under deforma- 
tions of r. In Ref. 1 (see also Appendix A) it is shown that the 
Euler-Lagrange equation for this variational principle implies 
that r is composed of pseudo-orbits tangential to the pseudo 
field 

= B - vV0xVf , (9) 

where v = n-B/n-VOxV^ is constant on each such pseudo- 
orbit. Intuitively, the pseudo field is constructed from the true 
field by adding a radial field that cancels the radial field caused 
by a perturbation away from a neighboring integrable field, 
while the poloidal field is unchanged [9]. Although QFMin 
curves in a Poincare section are not guaranteed to have the 
graph property, QFMin pseudo-orbits are easier to construct 
than ghost orbits and thus it would be advantageous to find a 
QFMin formulation whose pseudo-orbits are also ghost orbits. 

To construct a rational-rotational-transform QFMin surface 
we find periodic pseudo-orbits with rotational transform pjq. 
We vary v continuously over the range for which solutions can 
be found, the corresponding pseudo-orbits sweeping out rib- 
bons that may be joined to form the entire surface F. As the 
range includes v = 0, F includes the closed field lines, both 
action-minimizing and minimax, associated with the magnetic 
island chain with the given rotational transform. 

5. Comparison 

For illustration we use a model magnetic Hamiltonian (see 
Sec. 7), xo + X\7 consisting of an integrable part, xa - */'^/2, 



and a perturbation, 

X\ = ^Xm,n{^) cos(m0 - nO . (10) 

We use only two nonzero, (/^-independent, perturbation harmon- 
ics, X2,\ - 0.0010 and x^.i - -0.0005, to drive islands at the 
i = 1/2 and t - 2/3 rational surfaces. The degree of chaos 
induced by the perturbations is illustrated by a Poincare plot, 
shown with red dots in Fig. 2. Quadratic-flux-minimizing sur- 
faces (thick dashed lines) constructed as in Sec. 1 1 and La- 
grangian ghost surfaces (thin lines) constructed as in Sec. 12 as- 
sociated with 3 1 rationals between these two islands are shown. 
We term these QFMin surfaces uncorrected because no attempt 
at transforming the coordinate system to improve agreement be- 
tween QFMin and ghost surfaces has been made (i.e. we have 
taken - 9). On the scale of the figure, the quadratic-flux min- 
imizing surfaces and the ghost surfaces are indistinguishable. 

As the degree of chaos increases, the quadratic-flux min- 
imizing surfaces and the ghost surfaces appear to deviate, 
particularly those associated with the higher order rationals. 
Increasing the perturbation amplitude to X2,i - 0.0020 and 
^3,2 = -0.0010, the surfaces are shown in Fig. 3. 

The close agreement between QFMin curves and ghost 
curves in the case of moderate chaos for the coordinate choice 
used (action-angle coordinates in the unperturbed system) sug- 
gests a coordinate system in which they become identical may 
not be strongly perturbed away from action-angle coordinates. 
Below we point the way to achieving this unification. 

6. Modified Hamilton's Principle 

It has long been known [10, 2] that the quadratic flux (p2 
can be expressed in terms of action gradient, thus generalizing 
the QFMin concept to general Lagrangian/Hamiltonian systems 
and providing a link with the ghost surface approach. However, 
neither approach so far has provided a variational principle for 
individual pseudo-orbits. In this section we present a new and 
very useful variational principle for the QFMin pseudo-orbits, 
modifying Hamilton's Principle by adjoining the constraint 

0V^-dl= 0t/f = const. (11) 

c Jo 

This associates with a path r = r(f) an "area" ^ under the 
corresponding curve in the (^, 0) plane. By writing - 
{plq)C -1-00-1- 0(0' where is a periodic function averaging 
to and 0o is a constant approximately equal to the value of 
at which the pseudo-orbit cuts the section ^ = 0, we get 
J[ - 2nq&Q + 2n^pq. Thus, constraining, to be constant 
allows one to select a pseudo-orbit labelled by 0o, and 

r = rp,,/^|0o) , 0e[O,27r), (12) 

defines a family of paths covering a torus embedded in M^. 

To implement the constraint Eq. (11) using a Lagrange mul- 
tiplier V we define the pseudoaction 

Sy = S- vJ[, (13) 
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which is the same as the physical action Eq. (3) with A replaced 
by the vector pseudopotential Ay = A - y0V^. Thus, the 
verification of Hamilton's Principle for QFMin pseudo-orbits 
goes through exactly as that in Sec. 3 for physical magnetic 
field lines (except that we must use the restricted variations, 
6r = Ar, so that the endpoint contribution 2np6r-V^ vanishes): 
6Sy/6s = 6Sy/60 = implies rxBy = 0, where By is the 
pseudo magnetic field Eq. (9), as required for QFMin pseudo- 
orbits, for which r = By/B^ . 

7. Hamiltonian formulation 

By exploiting gauge freedom we may write 
A = (^V6'-x(i/',6',^)V^, for which S takes the famiUar 
form 



L 



{ijfde-xdo 



(14) 



This is the canonical form for the vector potential, as the equa- 
tions describing the field Unes are = dxldij/ and if/ = -dx/dO. 
These are Hamilton's equations, with the magnetic-field-line 
Hamiltonian. Henceforth we remove the arbitrariness in the ra- 
dial coordinate s by choosing it to be the flux function if/, so that 
s and 6 are now canonically conjugate phase-space variables for 
the given magnetic field, and the return map generated by field 
lines intersecting the Poincare section ^ = is area-preserving. 

From Hamilton's equations, we see that iix depends only on 
if/, e.g. X = Xoi>ff), then i// and 6 are action-angle coordinates 
[11]: (J/ is constant, and 6 increases linearly with f according 
to = *((/'). If the transform is rational, 6 = p/q for inte- 
ger p and q, there exists a continuous family of periodic field 
lines. Each periodic field line may be identified, for example, 
by the poloidal angle where the field line intersects the Poincare 
section ^ = 0. Together, the periodic field Unes form a ratio- 
nal surface. However, in a nonintegrable system no transfor- 
mation to action-angle coordinates exists and the Hamiltonian 
has the more general form;^; =xoiil/) +xi(^, 0,0- No matter 
how small xi is, the continuous family of periodic field lines 
is destroyed, being replaced by isolated action-minimizing and 
minimax orbits in an island chain. 

In canonical coordinates the QFMin pseudoaction, Eq. (13), 
is, from Eq. (14), Sy = §{\f/e-x - v0)Jf, up to 0(Ar), 



X/.nq 
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Setting the first variation to zero, we find Hamilton's equations 
for QFMin pseudo-orbits 



= d^Xv 
4' = -deXv . 



(16) 
(17) 



where ;ify =x + That is, v0 acts as a scalar pseudopotential 
that displaces pseudo-orbits poloidaUy, away from the periodic 
orbits found when y = 0. 



8. Lagrangian formulation 

An arbitrary trial curve, C, requires both the "position curve," 
6{^), and the "momentum curve," i/f(f), to be specified. 

It is advantageous to reduce the number of dependent vari- 
ables by transforming from the Hamiltonian phase-space de- 
scription to the Lagrangian configuration-space description in 
the standard way, eliminating the momentum ij/ in favor of the 
velocity given by Eq. (16). [We assume the "twist condition," 
X'q{'^) > 0, to allow unique inversion to give ij/ = 'J'y(d, 6, ^), so 
this transformation may not be possible for systems with non- 
monotonic rotational transform.] 

The velocity curve then defines the momentum curve, mak- 
ing Sy a functional of the position curve alone, while partially 
extremizing 



?y= r d^Ly{e,e,o, 
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where the pseudo-Lagrangian Ly = ^yO - XvQ^v, 0- (The 
physical field-line Lagrangian L and action S are obtained as 
the special case v = 0.) 

Now the total perturbation of the pseudoaction is of the form 
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where the Lagrangian pseudoaction gradient is 
6Sy^ ^dL^ _ d_ldLy\ 

66 86 d^\de}' 

9. Reconciling QFMin and Lagrangian Ghost surfaces 

We now seek to choose so that QFMin pseudo-orbits are 
also Lagrangian ghost pseudo-orbits as defined in Appendix B. 
First note that, to reconcile Eq. (16) and Eq. (47), we require 
= 0. That is, 



= @{9,0 



(21) 



Then *Py = Ly = L - v0, dLyld6 = dLjde, and members of 
our family of QFMin pseudo-orbits 



e = 0(f |0o) , 



(22) 



where 0o is an as yet arbitrary label [cf. Eq. (12)], satisfy the 
Euler-Lagrange equation 



6Sy 6S 

= — -y(0o)0,(0,f) = O, 
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(23) 



with @e{6,0 = d®(9,Old9. 

To reconcile QFMin and ghost orbits we require that the fam- 
ily of pseudo-orbits defined by Eqs (22) and (23) is the same 
family as is generated by Eq. (48). Thus the labels 0o and t 
must be functionally dependent: t = t(0()), c/t = T'(0o)<i0o- 
Eliminating 6S/66 between Eq. (23) and Eq. (48) and observ- 
ing that D&/D@o = iD6/D@o)@g we find the reconciliation 
condition 
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We now define ©o so that, for all ©o, 

t'(©o)v(0o) = -I, 

0(^|©o + 27r) = 0(^|©o) + 27r, 



where 6„ 



= 1 



deletes the resonant components 



(25) 



choosing //(^) so that Eq. (24) satisfies the solvability condition 
that the integral of both sides with respect to ©o over the interval 
[0, 2n] must be 2n, giving 
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10. Perturbative construction of QFMin-ghost surfaces 

For example, to approach this task perturbatively, consider 
the Lagrangian of the class corresponding to the case studied 
numerically in Sec. 5 



(which have been absorbed by vi — ^unlike the KAM problem, 

perturbation theory for almost-invariant tori is not inherently 
afflicted by small denominators, at least when * is a low-order 
rational). 

As the two Fourier coefficients in Eq. (31) are the same, 
0i(^|©o) = Qiitp^q^ + &o,0- However, ©i is not used in the 
calculation of di: to first order, ghost and uncorrected QFMin 
pseudo-orbits are identical, which, combined with the vanish- 
ing of m, is consistent with the near indistinguishability of 
these two almost-invariant surfaces in Fig. 2. 

At O(e^) we find 

I im'(m + m')^ V,„+ 
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L= ir-e/^ Vm,n expiimO - inO , 



(27) 



m,n=-oo 



with the reality condition = V-m-n (superscript •* de- 
noting complex conjugation), and e the expansion parameter. 
(The cases studied in Sec. 5 thus correspond to the choices 
Vi,! = ^^2,1/2, V3,2 = X3,2/2, and y„,„ = for [m,n} ^ 
{{2,1}, {-2,-1}, {3, 2}, {-3, -2}}.) 

As the unperturbed system is integrable, the expansions of v, 
fi, and 0(0, are of the form 



J" = 
= 



evi + e V2 + . . . , 
efii + e^H2 + ... , 



(28) 



where the prime on the sum over m' and n' indicates that the 
resonant terms, m'p = n'q, are to be deleted, and v^'„ = 
-iniVm,nSmp,nq- The term containing 0^'* gives the O(e^) dif- 
ference between ghost pseudo-orbits and uncorrected QFMin 
pseudo-orbits. 

We also see from Eq. (26) that //2 = -2€^{{d9\ld®of'), where 
(•) denotes averaging over ©o, so the correction factor jj. in the 
definition of ghost surfaces required for reconciliation does de- 
viate from unity by small amount, 0{e^). In the numerical work 
presented in Sec. 12 we take fi = I. 

11. Numerical construction of QFMin surfaces 



. . .) exp i{m0 - nQ , 



and of the {p, q) QFMin pseudo-orbits are of the form 
^(^I0o) = i;,.,f + 0o + 2«« + ^C + ---) 



Previously [9], periodic pseudo field lines were found by 
integrating the Hamiltonian equations of motion Eqs (16-17). 
The constrained action principle instead allows periodic pseudo 
orbits to be found variationally. This is of great benefit for 
strongly chaotic fields, as a defining characteristic of chaos is 
the exponential separation of initially nearby trajectories at a 
rate given by the Lyapunov exponent. Consequently, the intrin- 
X exp J - n^i^ + m©o] , sic numerical error associated with field-line-following meth- 

ods is guaranteed to increase as the trajectory becomes longer. 
^ / {m ipq-nj {^e6'^\ + e^6^^\ + . . The method of Lagrangian integration, also called variational 
m,n integration [12], avoids these problems. A Galerkin expansion 

of the trial curve. 



X exp i [(m tp^ - ^ -i- m©o] , (29) ' 



N-l 



where ip^q = p/q and i/' s dLjdO = 6 = 0f(f|0o). 
At 0(e) we find jUi = and 

oo 

Vi(0o) = -^imSmp,nqymflQ^Wm®o) , 



6= ir 



(33) 



(30) 



in, say, continuous, 2;r^-periodic basis functions m„ is in- 
serted into the action principle and the unknown ampUtudes a, 
are varied. 



where the Kronecker delta S„ipj,g selects Fourier coefficients 
such that mp = nq, resonant with (p, q) pseudo-orbits. 

The 0(e) contributions to the pseudo-orbit and pseudopoten- 
tial are given by 



66 ■■ 



N-l 

1=0 



(34) 



giving, from Eq. (19), 



J N-\ N-1 



5660' 



uj)aj + ... , 



(35) 



where the inner product (/, g), f{^) and g{^) arbitrary, is defined 
by 

{f,g)^ d^fg. (36) 
Jo 

Extremizing curves are then obtained, for example, by find- 
ing a zero of the A'^-dimensional gradient {6Sy/66,Ui) using 
Newton's method. This approach allows very high order pe- 
riodic orbits to be found, even for quite strongly chaotic fields 
[13]. 

For numerical work, we must describe O(^) with a finite set of 
parameters. It is simplest to use a piecewise hnear description, 
i.e. a finite-element expansion of 0(f) in tent functions: 

UiiO = i-[H(^-(i^i)H(^,-0(^-^.-i) 



(37) 



where H(-) is the Heaviside step function, f is to be evalu- 
ated mod 2nq to enforce periodicity [thus splitting the sup- 
port of Mo = between (0, fi) and (f/v-i. fw)]. fi = '^f. 
and Af = 2nq/N. At first, the piecewise-linear approxima- 
tion seems crude, but the great benefit is that the action integral 
can be calculated analytically in each interval [13]. The dis- 
cretized pseudoaction Sy becomes a rapidly computable func- 
tion of the independent parameters that describe the curve, 
[00,01, ... ,9i^-i], where 0, = ip/q)^i + a,. Extremal curves 
are found as zeros of the constrained action gradient, which 
are efliciently found using an A' -h 1 dimensional Newton's 
method. Note that the Hessian matrix {ui,{5^Syl6060)uj) as- 
sociated with the variations in the curve geometry is a cycUc 
tridiagonal matrix. 

12. Numerical construction of Lagrangian ghost surfaces 

We discretize Eq. (48) using the Galerkin method. That is, 
we substitute the ansatz Eq. (33) into Eq. (48) and project onto 
the finite basis {«,}, 



Dt 



6S 

~60 



(38) 



For a basis with global support the matrix (m,, Uj) is full and 
it must be inverted numerically to get Daj/Dr. For our finite- 
element basis the matrix is tridiagonal, 



{ui, Uj) - Af 



(Af)2 



(39) 



where the Kronecker is the identity matrix. The term 
(Si+ij - 25ij + 5,_ij)/(Af)^ is a finite-difference approximation 
to d^/di^^, so, assuming the action gradient is a smooth func- 
tion, the second term on the RHS is two orders in Af smaller 
than the first and so can be neglected. Thus Eq. (38) can be 
solved analytically to give 



Dt 



Dt 



1 



M,-, 



~60 



j_ds 

'A^d0i 



(40) 



13. Conclusion 

In this paper we have given generalized definitions of 
quadratic-flux-minimizing (QFMin) surfaces and ghost sur- 
faces, formulating them in terms of the magnetic field-line ac- 
tion. We have gone further by introducing a new constrained 
Hamilton's Principle for QFMin pseudo-orbits that a) facilitates 
a reconciliation between QFMin surfaces and ghost surfaces, 
and b) provides a better algorithm than o.d.e. integration for 
calculation of QFMin orbits in strongly chaotic regions. 

QFMin and ghost orbits have been computed for a model 
magnetic field perturbed away from an integrable system de- 
scribed in action-angle coordinates, and found almost, but not 
quite, to coincide. This is explained by finding a slight gen- 
eralization of our previous ghost surface definition that allows 
QFMin and ghost surfaces to be fully reconciled, in principle, 
the change in the ghost surfaces being O(e^). Using pertur- 
bation theory we find a choice of poloidal angle that unifies 
QFMin and ghost surfaces up to second order in e, with the 
difiference between ghost surfaces and uncorrected QFMin sur- 
faces being (9(e^). 

It remains for future work to implement this reconciliation at 
all levels of nonlinearity and to explore whether it is beneficial 
to use further generalizations of the ghost surface definition. 
Also requiring further investigation is the question of why ghost 
surfaces coincide so closely with temperature isosurfaces [4] 
for heat transport in chaotic magnetic fields. 

The innovations in this paper should also be applicable in 
general Ij-degree-of-freedom Hamiltonian systems and in the 
theory of area-preserving maps. 

A. QFMin Euler-Lagrange equation 

Here we rederive in canonical coordinates the Euler- 
Lagrange equation for quadratic-flux-minimizing surfaces 
found by Dewar et al. [1]. Consider a toroidal magnetic 
field in canonical form, B = Vx(4rV9 - x^O^ with Jaco- 
bian Ij ^ - Vt/^-V^xV^ . A toroidal surface is described by 
if/ = P(0, (). We define the tangential dynamics using the equa- 
tion - B'^/B^ as a constraint, and require the "radial" dynam- 
ics to be confined to the surface, so that tj/ - Pg0 + P^, where 
P„ = daP, a € {0,(\. The angle parametrization is arbitrary, 
and so let = ®{0,O a new poloidal angle. We define 
V as the projection, v s B • N, of the magnetic field onto the 
vector N = (e© -I- Pge^)x{<i( + P^e^) normal to T, where here 
(and only here) e©, are the derivatives of position with re- 
spect to 0, f at constant f, respectively: e© = dx/d@\( and 

= dx/d^\@. Combining these expressions, the tangential dy- 
namics may then be written = dx/difr and tfr - -dx/d0 - v@g, 
which may be recognized as the pseudo field, Eq. (9). The gen- 
eralized quadratic-flux functional Eq. (2) can be written 



^2 = ^JJflfedf (b-n)(b.n), 



(41) 



where N = (eg + Pge^)x(e( + P^e^), where eg = dx/d0\(, and 
now = dx/d^\g. The first variation in this functional due to 
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variations in the surface is 



C. Erratum 



5(p2 = J J ded^dP ^/g (B^de + d() v, (42) 

where we have used di ^lgB^) I da\^=p = d^(^^B'^)Pa + 
dai yfg^), for a e [9, f ) and (i e [0, to reflect flie constraint 
that when or ^ vary, must also vary to remain on the surface. 



The corrected version above includes the er- 
rata for Phys. Letts. A 373, 4409-4415 (2009), 
http ://dx.doi.org/l 0. 1 1 6/j .physleta.2009. 1 0.005 : 

On p. 4413, in the fifth line below Eq. (23), replace "t/r = 
T'{@o)de" with "dT = t' (00)6^00". 

On p. 4414, on the second line of Eq. (32), replace 



B. Generalized ghost surfaces 



In canonical coordinates Eq. (4) becomes, to 0((5r), 



ss 

Sr 



(43) 



XZTrq 
d^ [{e - d^x) 5<A - ('A + dex) SO + X^^] , 

where the second form follows directly from Eq. (14). Identify- 
ing coefficients of Sij/ and 66 between the two forms. 



6S 
6v 



(44) 



A sUghtly generalized definition of the Lagrangian definition 
[2] of ghost surfaces can be found by choosing the projection 
operator in Eq. (7) to be 



_ _ 1 1 

"ghost — ~^ilr^ilr + ) 



(45) 



where e — > is a switching factor required to obtain the La- 
grangian description, as described below, and //(^) = 0(1) is a 
factor we shall find necessary for unifying ghost and QFMin 
surfaces. Using Eq. (45) in Eq. (7), the generalized action- 
gradient flow deflning ghost surfaces in canonical coordinates 
becomes 



Dt 
DO 
Dt 



0-d^X 
e 

ip + dex 



(46) 



Next we take the e — > asymptotic limit, giving rise to two 
T-scales. On the short r-scale ^ adjusts exponentially fast to 
enforce 



(47) 



on the long r-scale. 

The action-gradient flow on the long r-scale defining ghost 
surfaces in Lagrangian form is thus 



DO _ _ I 6S 



(48) 



where 68166 = -tj/ - dgx = dL/dO - {d/dOidL/dO), with L = 
tj/O - X ™d ^ constrained to be a fimction of 6, 0, and f by 
Eq.(47). 



VI (0o) 



imtp,q - n) 



2 ' 



with 



(mtpg - nY 4-', 



and replace the sentence after this equation, 

"where the prime on the sum over m' and n' indicates that the 
resonant terms, m'p = n'q, are to be deleted." 

with 

"where the prime on the sum over m' and n' indicates that 
the resonant terms, m'p = n'q, are to be deleted, and v^), = 

iniVfii^fi6mp,nq • 
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